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^ ! Abstract 

m 

Environmental monitoring is often performed through a wireless sensor network, whose nodes are 
randomly deployed over the geographical region of interest. Sensors sample a physical phenomenon 
(the so-called field) and send their measurements to a sink node, which is in charge of reconstructing 
O . the field from such irregular samples. In this work, we focus on scenarios of practical interest 

where the sensor deployment is unfeasible in certain areas of the geographical region, e.g., due 
to terrain asperities, and the delivery of sensor measurements to the sink may fail due to fading 

ON : 

or to transmission collisions among sensors simultaneously accessing the wireless medium. Under 
\Q ■ 

these conditions, we carry out an asymptotic analysis and evaluate the quality of the estimation 
0^ ■ of a d-dimensional field (d > 1) when the sink uses linear filtering as a reconstruction technique. 

o ■ 

Specifically, given the matrix representing the sampling system, V, we derive both the moments 
and an expression of the limiting spectral distribution of VV H , as the size of V goes to infinity 
and its aspect ratio has a finite limit bounded away from zero. By using such asymptotic results, 
we approximate the mean square error on the estimated field through the 77-transform of VV H , and 
derive the sensor network performance under the conditions described above. 



I. Introduction 

Recently, a great deal of attention has been payed to wireless sensor networks whose nodes sample 
a physical phenomenon (hereinafter referred to as field), i.e., air temperature, light intensity, pollution 
levels or rain falls, and send their measurements to a central processing unit (or sink node). The sink 
is in charge of reconstructing the sensed field: if the field can be approximated as bandlimited in the 
time and space domain, then an estimate of the discrete spectrum can be obtained. 

However, the sensors measurements typically represent an irregular sampling of the field of interest, 
thus the sink operates based on a set of field samples that are not regularly spaced in the time and 
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space domain. The reasons for such an irregular sampling are multifold, (i) The sensors may be 
irregularly deployed in the geographical region of interest, either due to the adopted deployment 
procedure (e.g., sensors thrown out of an airplane ifflO. or due to the presence of terrain asperities 
and obstacles, (ii) The transmission of the measurements from the sensors to the central controller 
may fail due to bad channel propagation conditions (e.g., fading), or because collisions occur among 
the transmissions by sensors simultaneously attempting to access the channel. In this case, although 
the sample has been collected by the sensor, it will not be delivered to the central controller, (iii) The 
sensors may enter a low-power operational state (sleep mode), in order to save energy Q, 0. While 
in sleep mode, the nodes neither perform sensing operations nor transmit/receive any measurement, 
(iv) The sensors may be loosely synchronized, hence sense the field at different time instants. 

Clearly, sampling irregularities may result in a degradation of the reconstructed signal El. The 
work in Q investigates this issue in the context of sensor networks. Other interesting studies can 
be found in [|6j and 0, just to name a few, which address the perturbations of regular sampling in 
shift-invariant spaces and the reconstruction of irregularly sampled images in presence of measure 
noise (7J. 

In this work, our objective is to evaluate the performance of the field reconstruction when the 
coordinates in the (i-dimensional domain of the field samples, which reach the sink node, are randomly, 
independently distributed and the sensors measurements are noisy. We take as performance metric 
the mean square error (MSE) on the reconstructed field. As a reconstruction technique, we use linear 
filtering and we adopt the filter that minimizes the MSE (i.e., the LMMSE filter) QD-OID. The matrix 
representing the sampling system, in the following denoted by V, results to be a <i-fold Vandermonde 
matrixQ- By drawing on the results in O, ifTTI . we derive both the moments and an expression of the 
limiting spectral distribution (LSD) of VV H , as the size of V goes to infinity and its aspect ratio has 
a finite limit bounded away from zero. Then, by using such an asymptotic model, we approximate 
the MSE on the reconstructed field through the ^-transform |[T2l of VV H , and derive an expression 
for it. We apply our results to the study of network scenarios of practical interest, such as sensor 
sensor deployments with coverage holes, communication in presence of a fading channel, massively 
dense networks lfT3l . lfl4l . and networks using contention-based channel access techniques lfT5l . 

The rest of the paper is organized as follows. Section [II] reviews previous work, while Section UlTI 
describes the system model under study. In Section [IVl we first provide some useful definitions and 
introduce our performance metric, then we recall previous results on which we build our analysis. In 
Section IVl we derive asymptotic results concerning the moments and the LSD of VV H . Such results 

'An n x m matrix X is Vandermonde if its (i, j)— th entry, (X)y can be written as (X)y = Xj, i = 0, . . . , n — 1, j = 
1, ...,m. 
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are applied to different practical scenarios in Section [VI] Finally, Section IVIII concludes the paper. 

II. Related Work 

In the context of sensor networks, several works lfT6l - |fT9ll have studied the field reconstruction 
at the sink node in presence of spatial and temporal correlation among sensor measurements. In 
particular, in fl9l the observed field is a discrete vector of target positions and sensor observations 
are dependent. By modeling the sensor network as a channel encoder and exploiting some concepts 
from coding theory, the network capacity, defined as the maximum value of the ratio of the target 
positions to the number of sensors, is studied as a function of the noise, the sensing function and the 
sensor connectivity level. 

The paper by Dong and Tong EOl considers a dense sensor network where a MAC protocol is 
responsible to collect samples from network nodes. The work analyzes the impact of deterministic and 
random data collection strategies on the quality of field reconstruction. As a performance measure, 
the maximum of the reconstruction square error over the sensed field is employed, as opposed to our 
work where the mean square error is considered. Also, in EOll the field is a Gaussian random process 
and the sink always receives a sufficiently large number of samples so as to reconstruct the field with 
the required accuracy. 

The problem of reconstructing a bandlimited field from a set of irregular samples at unknown 
locations, instead, has been addressed in lf2TTl . There, the field is oversampled by irregularly spaced 
sensors; sensor positions are unknown but always equal to an integer multiple of the sampling interval. 
Different solution methods are proposed, and the conditions for which there exist multiple solutions 
or a unique solution are discussed. Differently from |2TI . we assume that the sink can either acquire 
or estimate the sensor locations and that the coordinates of the sampling points are randomly located 
over a finite <i-dimensional domain. 

As for previous results on Vandermonde matrices, in ifTTTl Ryan and Debbah considered a Van- 
dermonde matrix V with d = 1 and complex exponential entries, whose phases are i.i.d. with 
continuous distribution. Under such hypothesis, they obtained the important results that, given the 
phases distribution, the moments of VV H can be derived once the moments for the case with 
uniformly distributed phases are known. Also, a method for computing the moments of sums and 
products of Vandermonde matrices, for the non-folded case (i.e., d = 1), has recently appeared in ll22l : 
further insights on the extremal eigenvalues behavior, still for the case of non-folded Vandermonde 
matrices, can be found in |[23l . Moreover, in (9] it has been shown that the LSD of VV H converges 
to the Marcenko-Pastur distribution |[24l when V is d-fo\d Vandermonde with uniformly distributed 
phases and d — > oo. 

Note that, with respect to previous studies on Vandermonde matrices with entries that are randomly 
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distributed on the complex unit circle, in this work we obtain results on the LSD of VV where 
the entries of V have phases drawn from a generic continuous distribution. By relying on the results 
in (9], ifTTl . we show that such an LSD can be related to that of VV H when the phases of V are 
uniformly distributed on the complex unit circle. We also provide some numerical results that show 
the validity of our analysis. To our knowledge, these results have not been previously derived. We 
then apply them to the study of several practical scenarios in the context of sensor networks, although 
our findings can be useful for the study of other aspects of communications as well ifTTIl . 

III. Network Model 

We consider a network composed of m wireless sensors, which measure the value of a spatially- 
finite physical field defined over d dimensions, (d > 1). We denote by % = [— 5, the hypercube 
over which the sampling points fall, and we assume that the sampling points are i.i.d. randomly 
distributed variables, whose value is known to the sink node. Note that this is a fair assumption, 
as one can think of sensor nodes randomly deployed over the geographical region that has to be 
monitored, or, even in the case where the network topology is intended to have a regular structure, 
the actual node deployment may turn out to be random due to obstacles or terrain asperities. In 
addition, now and then the sensors may enter a low-operational mode (hence become inactive) in 
order to save energy, and they may be loosely synchronized. All the above conditions yield a set of 
randomly distributed samples of the field under observation, in both the time and the space domain Q. 

By truncating its Fourier series expansion, a physical field defined over d dimensions and with 
finite energy can be approximated in the region T~L as |[9l 

-^ /2 I>(^ TX (1) 

where n is the approximate one-sided bandwidth (per dimension) of the field, I = . . . , £d] T is 
a vector of integers, with = 0, . . . , n — 1, k = 1, . . . , d. The coefficient n~ d l 2 is a normalization 
factor and the function 

d 

maps uniquely the vector I over [0,n d -l]. a„(£) denotes the u(£)-th entry of the vector a of size n , 
which represents the approximated field spectrum, while the real vectors x g , q = 1, . . . , m represent 
the coordinates of the d-dimensional sampling points. In this work, we assume that x^, q = 1, . . . , m, 
are i.i.d. random vectors having a generic continuous distribution f x (z), z G %. In the specific case 
where x g are i.i.d with i.i.d. entries x q j,j = l,...,d, uniformly distributed in [—1/2, 1/2), we denote 
the distribution of x g by / u (z). 

The coordinates of the d-dimensional sampling points, however, are known to the sink node, 
because (i) either sensors are located at pre-defined positions or their position can be estimated 
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through a localization technique E51 . and (ii) the sampling time is either periodic or included in the 
information sent to the sink. 

Now, let s = [s(xi), . . . , s(x m )] T be the values of the samples at [xi,...,x m ], respectively. 
Following flU, El, we can write the vector s as a function of the field spectrum: 

s = /?n,™ /2 V H a (2) 

where V is the n d x m (i-fold Vandermonde matrix with entries 

(V„ ( *),,) = m' 1 ' 2 exp (-2vri€ T x g ) (3) 

randomly distributed on the complex circle of radius m -1 / 2 , and /3 n m is the ratio of the rows to the 
columns of V, i.e., 

B ~ n - 
Hn.m — 

m 

In general, the entries of a can be correlated with covariance matrix E[aa H ]. However, in the 
following, we restrict our attention to the class of fields characterized by E[aa H ] = <t 2 I. 

In the case where the sensor measurements, p = [pi, . . . ,p m ] , are noisy, then the relation between 
the sensor samples and the approximated field spectrum can be written as: 

p = s + n = /3-^ 2 V H a + n (4) 

where n is a m-size, zero-mean random vector representing the noise. Here, we assume a white noise, 
i.e., with covariance matrix E[nn H ] = <r 2 I m . Note that the additive white noise affecting the sensor 
measurements may be due to quantization, round-off errors or quality of the sensing device. 

IV. Preliminaries 

In this section, we report some definitions and previous results that are useful for our study. 
A. Useful definitions 

Let us consider an n x n non-negative definite random matrix A, whose eigenvalues are denoted 
by Aa,i, • • • , AA,n- 

Definition 4.1: The average empirical cumulative distribution of the eigenvalues of A is defined as 
^\i( z ) = n 2~27=i ^ [H^A.i ^ z }]> where the superscript (n) indicates that we refer to a system with 
size n and 1{-} is the indicator function. If (z) converges as n — > oo, then lim n _ s . 0O (z) = 
F\ A (z). The corresponding limiting probability density function, or limiting spectral distribution 
(LSD), is denoted by /a a (-)- 

Definition 4.2: The ^-transform of A is given by: 



??A l) (7)=E[tr{( 7 A + I)- 1 } 



E 



1 n 1 



n 



(5) 



September 6, 2010 



DRAFT 



where tr{-} is the normalized matrix trace operator and 7 is a non-negative real number. If r]£ (7) 
converges as n — > 00, then the corresponding limit is T]a(i) = E[(7Aa + 1) ] lfl2l p. 40], where Aa 
is the generic asymptotic eigenvalue of A, whose distribution is f\ A (z), and the average is computed 
with respect to A a Ifl2ll . 

Next, consider the matrix V as defined in (O and that the LMMSE filter is used for field 
reconstruction. Then, the estimate of the unknown vector a in (01), given y and V, is obtained 
by computing a = E[ap H ]E[pp H ]~ 1 p. Through easy computations and using the Sherman-Morrison- 
Woodbury identity, we can obtain the MSE as 

MSE (n) = <7~ 2 E 



tr|(^ 2 /3^VV H +a- 2 l) '} 



\Pn,m / 

where 7 = cr^/cr 2 denotes the signal-to-noise ratio on the sensor measurements, and we employed 
the definition of the 77-transform given in (0. 

Next, we approximate the MSE of the finite size system in (0]) through an asymptotic model, which 
assumes the size of V to grow to infinity while the ratio of its number of rows to its number of 
columns tends to a finite limit, j3, greater than zero, i.e., we assume 

lim P nm = P 

n,m—> 00 

Indeed, in our recent works ll8Tl — ifTOTl it was shown that this asymptotic model provides a tight 
approximation of the MSE of the finite size system, already for small values of n and m. Under 
these conditions, we therefore define the asymptotic expression of the MSE as ifTOl : 

MSEoo = lim MSE (n) = WV " {l/P) (V) 

n,m— >■ 00 

if the limit exists. 



B. Previous results 

Vandermonde matrices have been studied in a number of recent works lf8ll — iTTTTl - Specifically, (9] 
considered the case where the vectors x g are i.i.d., for q = 1, . . . , m, and their entries, x g j are i.i.d. 
random variables with uniform distribution in [—1/2,1/2). The work there studied the eigenvalue 
distribution of VV H for both finite and infinite (i.e., m, n — > 00) matrix size. Although an explicit 
expression of such LSD is still unknown, Q provided an algorithm to compute its moments of any 
order in closed form. 

Indeed, as n, m — > 00 with j3 n ^ m = n d /m having a finite limit j3 > 0, in (9] it was shown that the 
p-th moment of the generic asymptotic eigenvalue of VV H , denoted by A, is given by 

k=i u>efi P ,k 
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where f\^ u (d,f3,z) represents the distribution of A. Moreover, $7 p & is the set of partitions of the 
set V = {1,2, ... ,p} in k subsets, and v(cj) G (0, 1], u £ ^l p ^ is a rational number that can be 
analytically computed from u> following the procedure described in (9l- The subscript u in M p ^p, u 
and f\^ u (d, ft, z) indicates that a uniform distribution of the entries of x,j is considered in the matrix 
V. 

In |9} it was also shown that when n,m,d — > oo, with /3 n m = n d /m having a finite limit (3 > 0, the 
eigenvalue distribution f\^ u (d, (3, z) converges to the Marccenko-Pastur law Il24l . A similar result ifTOll 
also applies when the vectors x g (q = 1, . . . ,m) are independent but not i.i.d., with equally spaced 
averages. 

More recently, Ryan and Debbah in flTTl considered d = 1 and the case where the random variables 
x q i, q = 1, . . . , m, are i.i.d. with continuous distribution f x (z), < z < 1. Under such hypothesis, 

it was shown that the asymptotic moments of VV H can be written as 

v 

M pXP , x = Y^hP p ' k £ v(u) (8) 
k=i a;ef2 p , fc 

where the terms 1^ depend on the phase distribution f x (z) and are given by 

Jo 

for k > 1. The subscript x in M Pi i j( g :X indicates that in the matrix V the random variables x q ± have a 
generic continuous distribution f x (z). Note that for the uniform distribution we have Ik = 1, for all 
k. The important result in ([8]) states that, given f3, if the moments of VV H are known for uniformly 
distributed phases, they can be readily obtained for any continuous phase distribution f x (z). 

V. Vandermonde matrices with generic phase distribution 

In this work, we extend the above results by considering a sampling system defined over d > 1 
dimensions with nonuniform sample distribution, where samples may be irregularly spaced in the 
time and spatial domains, as it occurs in wireless sensor networks. Being our goal the estimation of 
the quality of the reconstructed field, we aim at deriving the asymptotic MSE (i.e., ??vv H (7//3))- 

We start by considering a generic continuous distribution, f x (z), z G T~L of the samples measured by 
the sensors over the d-dimensional domain. We state the theorem below, which gives the asymptotic 
expression of the generic moment of VV H , for d > 1. 

Theorem 5.1: Let V a (f-fold n d x m Vandermonde matrix with entries given by (0) where the 
vectors x g , q = l,...,m, are i.i.d. and have continuous distribution f x (z). Then, for n, m — > oo, 
with fin^m = n d /m having a finite limit j3 > 0, the p-th moment of VV H is given by 

M Mx = jrp p ~ kl k E v ^ d (9) 
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where = f n f x (z) k dz and the terms v{(v) are defined as in 19] . 
The proof is given in Appendix |A) 

Using Theorem [5d1 and the definition of Ik, it it possible to show the theorem below, which provides 
the LSD of VV H . 
Theorem 5.2: Let 

• V be a d-fold n d x m Vandermonde matrix with entries given by ((3]) where the vectors x 9 , 
q = 1, . . . , m, are i.i.d. and have continuous distribution f x (z), z £ Ti 

• A be the set where f x (z) is strictly positive, i.e., A = {z G H\f x (z) > 0} 

• the cumulative density function 

G x (y) = j^\{zeA\f x (z)<y}\ (10) 

defined^ for y > and let g x (y) be its corresponding probability density function. 
Then, the LSD of VV H , for n, m — > oo with f3 n , m = n d /m having a finite limit /3 > 0, is given by 

/A,x(d,^2) = (l-|^|)^) + W r —fx,uUAA dy (11) 

Jo y \ y yj 

Proof: The proof can be found in Appendix |B] ■ 
From Theorem 15.21 the corollary below follows. 

Corollary 5.1: Consider f x (z) such that f x (z) > Vz G Then, let us denote by / X '(z) a scaled 
version of this function, so that 

/,(«) = ( (12) 
[ zeU\C 

where C C Ti. It can be shown that 

/a^M = (l-c),5(z) + C 2 /A,x(d,c/3, C z) (13) 

where c = |C|. 

Proof: The proof can be found in Appendix |C] ■ 
As an example of the result given in Corollary 15.11 consider that a unidimensional (d = 1) sensor 
network monitors the segment % = [—1/2,1/2]. Due to terrain irregularities and obstacles, nodes 
are deployed with uniform distribution only in the range [— c/2, c/2] (with c G [0, 1)). We therefore 
have f X '(z) = 1/c for —c/2 < z < c/2 and elsewhere. Moreover, f\^ x {l, (3, z) = f\ >u (l,l3, z). 
The expression of f\ )X >(l, (3, z) is given by (fL3l . by replacing d = 1 and the subscript x with the 
subscript u. 

This result is well supported by simulations as shown in Figures |l(a) and |l(b)| In the plots, 



(n) In) 

we compare the asymptotic empirical spectral distribution (AESD) f^ , and /I J instead of the 

2 IS denotes the measure of the set 5 
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(a) 



(b) 



Fig. 1. Comparison between the curve representing fj™ ,{l,/3n,m,z) and the one representing the empirical function 

c 2 fi n i(d,P„, m c,zc). (a) n>m = 0.8, c = 0.8 and (b) n>m = 0.2, c = 0.5. 



LSDs fx )X > and /x jU since an analytic expression of f\ jU is still unknown. However, in ll8l- lfT0l it is 

(n) 

shown that, already for small values of n, the AESD k ' appears to rapidly converge to a limiting 
distribution. Figure |l(a)| refers to the case (3 n ,m = 0.8 and c = 0.8. The solid and dashed lines 
represent, respectively, the functions /j^,(l,/3, z) and c 2 /^(l, (3 niTn c, zc), for n = 100. Note that 
the probability mass of /^(l, /3 n m , z) at z = is not shown for simplicity. Similarly, Figure [1(5)1 
shows the case /3 n m = 0.2 and c = 0.5. As evident from these plots, the match between the two 
functions is excellent for any parameter setting, thus supporting our findings. 

Since we are interested in evaluating the MSE, taking into account the result in (0), we now apply 
the definition of the ^-transform to (fTTT) . The corollary below immediately follows. 

Corollary 5.2: The r/-transform of VV H is given by 

»k(d,0,7) = l~\A\ + \A\ jf g x (y) nu (d, ^, 7 y) dy (14) 
hence, the asymptotic MSE on the reconstructed field, defined in (|7]), is given by 

MSE^ = Vx (d, p,l) = l- \A\ + \A\ g x (y)riu U ^ dy (15) 

Proof: The proof can be found in Appendix |E] ■ 
In (fT4l . in order to avoid a heavy notation we referred to ?7vv H (^> Pil) as "Ux{d, /5, 7) when the 
phases of the entries of V follow a generic random continuous distribution, while r] u (d, {3, 7) refers 
to the case where the phases are uniformly distributed. 

Remark 5.1: Since g x (y) > and r/ u (d, f3/y, jy/fi) > 0, the integral in the right hand side of (fT4l 
is positive, then rj x (d, ^,7/ /3) > 1 — |«4|. It follows that the MSE is lower-bounded by the measure 
of the total area where the probability of finding a sensor is zero. This clearly suggests that, in order 
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to obtain a good quality of the field reconstructed at the sink node, this area must be a small fraction 
of the region under observation. 

Next, we observe that, in the case of massively dense networks where the number of sampling 
sensors is much larger than the number of harmonics considered in the approximated field, i.e., /3<1, 
an interesting result holds: 

Corollary 5.3: Let A be the set where f x (z) is strictly positive; then 

lim f XjX (d, (3, z) = {l-\A\)5(z) + \A\g x {z) (16) 

Proof: The proof can be found in Appendix [D) ■ 
Thus, as evident from Corollary 15.31 for the limit of /3 — > 0, the LSD of VV H is the density of the 
density of the phase distribution /^(z). 

Furthermore, for massively dense networks, we have: 

Corollary 5.4: Let A be the set where f x {z.) is strictly positive; then 

Hm^d, 0,7/0) = l-l-^l (17) 

Proof: The proof can be found in Appendix 10 ■ 
Remark 5.2: The result in (fTTT ) shows that even for massively dense networks 1 — \A\ is the 
minimum achievable MSEqo, when an area A cannot covered by sensors. 

VI. Field reconstruction in presence of losses 

Here, we provide examples of how our results can be used in wireless sensor networks to investigate 
the impact of a random distribution of the coordinates of the sampling points on the quality of the 
reconstructed field. In particular, we first consider a wireless channel affected by fading, and then the 
effects of contention-based channel access. 

A. Sensor network performance with fading communication channel 

We consider a wireless sensor network whose nodes are uniformly distributed over a geograph- 
ical region. Without loss of generality, we assume a square region of unitary side (d = 2, % = 
[— 1/2, +1/2] 2 ), where the sink is located at the center and has coordinates (z±, z%) = (0, 0). Through 
direct transmissions, the sensors periodically send messages to the sink, including their measurements. 
At every sample period, a sensor message is correctly received at the sink if its signal-to-noise ratio 
(SNR) exceeds a threshold r. The communication channel is assumed to be affected by slow fading 
and to be stationary over the message duration. 

Let d be the distance between a generic sensor and the sink. Then, the signal to noise ratio at the 
receiver is given by 

SNR(d) = s\h\ 2 d~ 2 

September 6, 2010 DRAFT 
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where h ~ A/c(0, 1) is a circularly symmetric Gaussian complex random variable representing the 
channel gain, and s is the signal to noise ratio in the absence of fading and when the sensor-sink 
distance is d = 1. 

The probability that a message is correctly received at the sink is given by 

P(SNR(d)>r) = P i\h\ 2 >T— 



\r-) 



s J 

exp (-ad 2 ) (18) 



2 



with a = t/s and F\^(z) = 1 — e~ z being the cumulative density function of \h\ 

The probability density f x (z\,z 2 ) corresponding to sensors at distance d = \J z\ + z\, — 1/2 < 
zi, z 2 < 1/2 from the sink and successfully sending a message is then given by 

f u ( Zl , z 2 )F (SNR (^zf+zfj > r) 

fx&l, z 2) = 7 7 ; x 

// w /„(^i I ^)P(SNR[v^ + ^) > t ) dz i dz 2 
where f u {z\, z-z) = 1 Vzi, z 2 , is the density representing the sensor deployment (recall that nodes are 
assumed to be uniformly distributed in the region hence their density is constant and equal to 1). 
Using (fT8l) . we obtain: 

exp (-a(zl + zl)) 



fx(zi,Z 2 ) 



ff n exp (-a(zf + zl)) dz\ dz 2 

b^ V (-a(zl + z 2 2 )) (19) 



where 

b 



1 = J J exp (— a(zf + z 2 )) dz\ dz 2 = —erf 2 

In order to compute (fT4l) . we need the function g x (y), i-e-, the density of f x {z\,z 2 ). Note that 
fx(z±, z 2 ) is circularly symmetric with respect to (z±, z 2 ) = (0, 0). Let y be the value of density of the 



sampling points at distance d y = -y/ ' z\ + i| from the sink. Then, from ( fT9l we obtain d y = y ^ log |, 
thus the network area where the density is lower than y is given by 

G x (y) = 1 - ndl 

for < d y < 1/2, i.e., 6e~ a/4 < y < 6. For 1/2 < d y < y/2/2, it is possible to show that 

G x (y) = 1 - ^/Ad 2 - 1 - 4 (tt - 4 cos" 1 



In conclusion, 



G x {y) 



1 y/6 > 1 

l-ird^ e~ a / A <y/b<l 

1 - ^Ju 2 y -l-d 2 y (tt - 4 cos" 1 g^) e~ a / 2 <y/b< e~ a / 4 

y/b < e~ a / 2 



(20) 
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and 



9x{y) 



• in 



y/b>i 

e -«/4 <y/b<l 



i^-^os- 1 ^) e-/ 2 < y/b < 
y/6 < e" a / 2 

Since |^4| = 1, then the asymptotic MSE can be obtained by computing 

MSE^ = Vx (2, 7 //3) = f g x (y)r, u {% £ ^ dy 

Je-°/ 2 V y P J 



(21) 



(22) 
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Fig. 2. Transmissions in presence of fading: density g x {y) for a = 0,5,10 dB, i.e., for different values of the SNR 
threshold r. 

Figure [2] shows the density g x (y) for a = 0,5,10 dB. Note that a = t/s, thus for a fixed s (i.e., 
the signal to noise ratio at distance D in the absence of fading) the parameter a is proportional to the 
SNR threshold r. In particular, as r decreases, the probability that a message successfully reaches the 
destination increases and, thus, the spatial distribution of correctly received samples, f x (zi, z?), tends 
to the uniform distribution f u {zi,Z2). As a consequence, the density of fx{z\,Z2), i.e., g x {y), for 
a = and 5 dB is concentrated close to y = 1. However, for high values of r, messages originated 
from sensor nodes located far from the sink are successfully received with low probability. Thus, 
g x (y) shows a significant probability mass around y = 0. 

Figure [3] shows the effect of the fading channel on the MSE of the reconstructed field (dashed 
lines), and compares the obtained results with the MSE obtained in absence of fading (solid lines). 
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Y(dB) 

Fig. 3. MSE of the reconstructed field in absence (f u (z)) and in presence (f x (z)) of fading, as the signal to noise ratio 
on the sensor measurements varies. 

The plot considers different values of (3, namely, /3 = 0.2,0.4,0.6,0.8, and a = 5 dB. The MSE 
is plotted versus the signal to noise ratio on the sensor measurements, 7. The curves have been 
obtained by numerically computing (1221 ). where g x (y) is given by (|2TT > and rj u (2, f3,j//3) is replaced 
by r]^ (2, /3 n>m , 7//3 n ,m)> wrtn n = 10. Recall that the analytic expression of the LSD f\, u (d, (3, z) is 
unknown, hence in the numerical results we considered the AESD f^(d, /3 n , m , z) instead. We observe 
that for low values /3, in spite of the presence of fading, the sink node still receives a large number of 
samples from the sensors, hence the degradation of the MSE shown in Figure [3] is negligible. On the 
contrary, for /3 > 0.4 (i.e., for a larger value of the ratio of the number of harmonics composing the 
approximated field to the number of sensors), the reconstruction performance degrades significantly 
and this is particularly evident in presence of high values of 7. 

In the case of massively dense networks, the LSD of VV H is given by (fl6l ) and from (fTTT ) we know 
that the MSE tends to as (3 — > 00. This result is confirmed by the plot in Figure |4] which shows 
the AESD f^(2, (3 n ,m, z), for \A\ = 1, a = h dB, and n = 10. The behavior of such a function is 
compared with the density g x (z) as (3 varies. We note that, as (3 decreases, the matching between 
/jf2(2, P n ,mi z ) and g x (z) improves, and the latter represents an excellent approximation already for 
f3 = 0.01, as predicted by the result in (15.31 ). 
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Fig. 4. Massively dense networks: empirical function /} n j(2,/3,z) in presence of fading, with a = 5 dB and \A\ = 1. 
The curves obtained for different values of /3 are compared with the density g x (z). 



B. Measurements gathering through contention-based channel access 

In environmental monitoring applications, it is often desirable to vary the resolution level with 
which the field measurements are taken over the region under observation, depending on the field 
variations and the interest level of the different locations Q, 0. It follows that the number of samples 
generated by the sensors network (i.e., the offered traffic load) varies in the spatial domain. 

To represent such a scenario, we consider a wireless sensor network whose nodes are uniformly 
deployed over a square region. We also identify L areas, Aj i = 1, ... ,L, each corresponding to 
a different value of the offered traffic load. As often assumed in the literature (see e.g., |[26ll . E710 
and widely applied in the practice, the network is divided into clusters and a hierarchy of clusters is 
created. More specifically, at the first hierarchical layer, layer 1, the sensors are grouped into clusters, 
each of which is controlled by a cluster-head. The cluster-head is in charge of handling all traffic 
packets it receives from the nodes. At a given layer h > 1 of the hierarchy, the cluster-heads are 
grouped into clusters on their turn and forward the traffic to their parent cluster-head. At the highest 
layer, layer H, we have only one cluster whose cluster-head coincides with the sink node. Without 
loss of generality, we assume that the cluster at the H-th layer is composed of L cluster-heads, each 
handling the traffic generated within one of the L areas defined above. 

As for the medium access control (MAC) layer, we consider that the nodes implement the IEEE 
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802.15.4 standard specifications for wireless sensor networks |[T5l . In particular, all nodes within a 
cluster are in radio visibility of each other and use the slotted carrier-sense multiple-access/collision 
avoidance (CSMA/CA) technique lfT31 . This is a contention-based scheme and transmissions may 
fail if two or more sensors access the channel at the same time. Inter-cluster interference is instead 
avoided by assigning different frequency channels to neighboring clusters. We consider that packets, 
whose transmission fails, are discarded. 

In order to derive the probability that a packet transmission fails within a cluster due to collision, 
we use the Markov chain model presented in |[28l . We denote by m^h the average number of sensors 
belonging to the generic cluster at the h-th layer of the hierarchical architecture, in area Ai (i = 
1, . . . , L and h = 1, . . . , H). Similarly, we define Aj^ as the average traffic load per node, again 
within the generic cluster at the h-th layer, in area Ai. Then, we set the size of the packet pay load to 
32 bytes, and the value of the other parameters as in |[28l . Under this setting, we compute the value 
of the collision probability within the generic cluster at layer h, in area Ai, as a function of m^h 
and X^h, i.e., P c (i,h) j28l . Furthermore, we observe that at the generic layer h, with 1 < h < H, 
a node, which acts as cluster-head at layer h — 1 in area A4, will have a traffic load equal to 

\,h = mi,h-l\,h-i [1 - PS, h - 1)]. 

It follows that the probability that a packet is successfully delivered to the corresponding /i-layer 
cluster-head within area A4 (i = l,...,L) can be obtained as P s (i,h) = 1 — P c (i,h). Then, the 
probability that a measurement generated by a sensor located in A4 (i = 1, . . . , L) is successfully 
delivered to the sink is given by: 



Next, denoting by \A-i\ the measure of A4, we define 



as the normalized probability that a message is successfully delivered to the sink. Then, the spatial 
density of the sensors successfully sending their message is as follows: 

fx(zi,z 2 ) = p s (i) V(zi, z 2 ) £ Ai, i = l,...,L. 

The density of f x (z±,Z2) is therefore given by 

L 

9x{y) = ^2 \ Aj -\ 5 (v-pS)) ■ 

i=l 

and the asymptotic MSE is given by 



H 



P s (i) = l- \P B (i,h). 



h=l 




i=l 
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Fig. 5. Comparison between the case where transmission collisions are taken into account (f x (z)) and the case where 
all measurements successfully reach the sink (/„(z)). The MSE is shown as a function of f3 and for different values of 
signal-to-noise ratio (H = 3, \Ai\ = 1/4, Vi, Ai,i = 1(T 3 , A 2 ,i = 2 ■ 1(T 4 , A 3 ,i = 2 • 1(T 4 , A 4 ,i = 2 ■ 10 _s ). 



Figures [5] and [6] show the impact of collisions due to the contention-based channel access, on 
the quality of the reconstructed field. In particular, they compare the MSE of the reconstructed field 
when collisions are taken into account (/ x (z)) with the one obtained in the idealistic case where all 
messages (measurements) sent by the sensors successfully reach the sink (f u (z)). The results refer to 
a square region of unitary side, where there are four areas of equal size (\Ai\ = 1/4, i = 1, . . . ,4) but 
corresponding to different resolution levels in the measurements collection (i.e., they are characterized 
by different traffic loads); the number of hierarchical levels is set to H = 3. We set Ai,i = 10 -3 , 
A 2 ,i = 2 • 1(T 4 , A 3 ,i = 2 • 1CT 4 , A 4) i = 2 • 1(T 5 in Figure and a higher traffic load in Figure H 
i.e., Ai,i = 5 • lCT 3 , A 2 ,i = 1CT 3 , A 3 ,i = 1CT 3 , A 4 ,i = 10" 4 . 

Looking at the plots, we observe that both j3 and 7 have a significant impact of the obtained MSE, 
with the MSE increasing as /3 grows and smaller values of 7 are considered. Most interestingly, by 
comparing the two figures, we can see that as the traffic load, hence the collision probability, increases, 
the performance derived taking into account the contention-based channel access significantly differs 
from the idealistic one. Furthermore, the latter effect is particularly evident as 7 increases, since the 
higher the signal-to-noise ratio, the more valuable the samples sent by the sensors toward the sink. 
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Fig. 6. Comparison between the case where transmission collisions are taken into account (f x (z)) and the case where all 
measurement transmissions are successful (/«(z)). The MSE is shown as /3 varies and for different values of signal-to-noise 
ratio (H = 3, \At\ = 1/4, Vi, A M = 5 • lO -3 , A 2 ,i = 1(T 3 , A 3 ,i = 1(T 3 , A 4 ,i = 1(T 4 ). 



VII. Conclusion 

We studied the performance of a wireless network whose nodes sense a multi-dimensional field 
and transfer their measurements to a sink node. As often happens in practical cases, we assumed 
the sensors to be randomly deployed over (the whole or only a portion of) the region of interest, 
and that their measurements may be lost due to fading or transmission collisions over the wireless 
channel. We modeled the sampling system through a multi-folded Vandermonde matrix V and, by 
using asymptotic analysis, we approximated the MSE of the field, which the sink node reconstructs 
from the received sensor measurements with the r/-transform of VV H . 

Our results clearly indicate that the percentage of region where sensors cannot be deployed must 
be extremely small if an accurate field estimation has to be obtained. Also, the effect of losses due to 
fading or transmission collisions can be greatly mitigated provided that a suitable value for the ratio 
between the number of harmonics approximating the field bandwidth and the number of sensors is 
selected. 
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Appendix A 
Proof of Theorem [57T] 

The p-th moment of the asymptotic eigenvalue distribution of VV H can be expressed as [9[ 



Mp,d,p, x = lim E 

n,m— >oo 



tr 



{(w-n 



where tr{-} is the normalized matrix trace operator. The matrix power can be expanded as a multiple 
sum over the entries of V: 

1 



Mp,d,fs, x = lim 

n,m— >oo mPn 



V V IE e j27r£lT(x 'f~ x ' ! i ) ---e" j27r£pT(x ' J p- 1_x * ) 
?i,...,£ p qi,...,q P 



where gi, • • • , g p , gj = 1, . . . , m are integer indices and £\, . . . , i p , t- L = [£i t i, . . . , ti^] 7 , (4,3 
0, . . . , n — 1 are the indices identifying the rows of V. Since, 



n-l 



+ -+ii,dX d ) 



TT- 

11 l _ e j 271-2^ 



for i = 1, . . . ,p and the elements of x are i.i.d., we have that 

1 



M 



p,d,(3,x 



lim 

n,m— >■ 00 



mPn d 



E E 



1 7=1 1 - 



\2n(x a . i—x 



1=1 3 



where the index i is to be considered modulo p, i.e., p + 1 = 1. As for the sum over the indices 
gi, . . . , q p we note that any choice of q = [qi, ■ ■ ■ ,q p ] T induces a partition u> of the set, T 3 = {1, ■ ■ ■ ,p} 
in k subsets Vi, ■ ■ ■ , Vk, 1 < k < p, under the equality relation [9]. In the following, we denote by 
Op fc the set of partitions of V in k subsets, 1 < k < p. Since there are m k possible vectors q inducing 
a given partition u G 0„ we can write the p-th moment as 



M, 



lim y y 

n,m— >oo — J ^— 4 



np-k 
Hn.m 



k=i u>efi Ptk 
p 



n d(p—k+l) 



E 



V d 1 \2Trn(xu 

nn 



3 -^^i+i ,3 , 



1.7=1 



!™ E^r™ E 

n,m— >oo ^— » ' z — J 



E[* w (xi J ...,x fc )] 



fc=i uiefip.k 



11 



d(p-k+l) 



(23) 



where $ w (x x , . . . ,x fc ) = n|=i -^(^y, ■ ■ ■ ,^fej), F w (xij, . . . , = nf=i 1 6 j27r / "" J _ x . 
and cjj S {1, . . . , k} is the index of the subset of V containing i. Recall that p + 1 = 1 and that 
Pnm = n d /m. Moreover, since the vectors x are i.i.d., we removed the dependence on the subscript 



g- 



Following the same steps as in |[TT1 Appendix H], we compute the limit 



E[$ w (xi,... ) x fc )] 



lim 

n—¥oo fi 



d{p-k+l) 



lim / /x(xi)---/a.(Xfc 



$ w (xi, • • • ,Xjfc) 



<f(p+l-fc) 



dxi ■ ■ ■ dx fc 



lim / /x(xi)---/ x (x fc )T 
n— >-oo Iqik J - J - 

J=l 



i^o;(^lj) • • • j •Ekj) 



+l-k 



dxij • • • dx 



kj 
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We then define = [xhi, Yh] where = [x^ • • • > z/id]> for /i = l,...,k and we integrate first 
with respect to the variables xn, . . . , x k \ obtaining 



E[$ w (xi,...,x fc )] v f n i .yj F u (xij,...,x kj ) 

n— too n a(p—K+i) n— >oo J\_i iiw-ijfe - LJ - n" + 

I 2 ' 2 J j = 2 

where 

Gu(yi, ■ ■ ■ ,Yk) = [ F " Xn ^ k ll Xkl fx([xn, yi]) • • • /x([gfci) yfc]) dsn • • • dxfci (24) 
In HU Appendix H] it was shown that, because of the properties of F u (xu, . . . ,Xk\), 

lim f 11 ' , ' — — dxu • • • dx/d = 

where 

B e = {(xn, . . . ,x fc i)||x w - xa| > e,for some h,£} 



for any e > 0. This means that the integral in (1241 ) can be limited to the x\\, . . . , x k \ on the diagonal 
where in = • • • = x k \. Therefore 

lim G w (yi,...,y fc ) = lim / F " %kl f x ([x kl , yi]) • • • f x ([x k i, Yk]) dsn ■ ■ • dx fc i 

ra— too n— yoo Jy_i ljfe 71^ ^ 

L~2'2J/l=l Vl _ 3'3J h = l / 

= / T /a;([xfci,y^)dxjfci (25) 



Note that the limit 



»(«) = iim /" F - (xil 'v;; Xfcl) T7 dx w 

J t I'll h=l 
does not depend on x k \ and the coefficient v(u>) £ [0, 1] is described in 

Next, iterating this procedure by integrating over the variables, xy, . . . , xjy, j = 2, . . . , d we finally 

get 

) d / fx(xki, x k d) k dxfci • • • dx kd 



lim ^ = 



2 ' 2 J 



v(u) d [ /,(x fc ) fc dx A 
JH 



v(co) d I k (26) 



where we defined I k = J n f x {~x-k) k dx^. It follows that 



k=i u>efi Ptk 

which proves the theorem. Note that when the entries of x 9 = [x 9 i, . . . ,x q d] T are independent with 
continuous distribution f x j(zj) such that / x (x) = Ylj=i fx, ( x j)> we have I k = Y\.j=i ^k,j with 
= J[-i/2,l/2] $ x A x ) k d^- 
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Appendix B 
Proof of Theorem [572] 

From Theorem 15. II and the definition of Ik, we have that 

m Mx = j^^ k h e <") d 

k=l wESlp.fc 
V 



f i2^- k f x (z) k Y v(") d dz. (27) 

U k=l u6fi„ t 



Next, we define the set A where f x (z) is strictly positive as 

A = {z£H\f x (z)>0} 
Note that for z £ % \ A the contribution to the integral in (l2"7T ) is zero. Thus, 

/ j^^f^t E v ^ dz 



k=l wEO- 



fc=i ueS! P , t 

fc=i wen„,t 
= / /^zfM^^dz (28) 

where for any z G A, M p ^>u^ is the p-th moment of VV H when the phases are uniformly distributed 
in % and the ratio ft' iz) is given by 



Note also that (|28T ) holds for p > 1 since, by definition, the zero-th moment of any distribution is 
equal to 1. Expression (1281 ) allows us to write the moments of VV H for any distribution f x (z), given 
the moments for uniformly distributed phases. Likewise, it is possible to describe the LSD of VV H , 
for any continuous f x (z), in terms of the LSD obtained for uniformly distributed phases. Indeed, let 
us denote the Laplace transform of f\ >x (d, ft, z) by L\ x (d, ft, s) if it exists. Then, whenever the sum 
converges 

°° s p 

L x , x {d, ft, s) = E —M Pt d,p,x 
p=0 p ' 
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Since Mo d^g x = 1 for any distribution, using (1281 ) we obtain: 

L XjX (d,p,s) = 1 + Yl^i J A fx(z) P M pA p, (zlu dz 

= l-\A\ + \A\+ J^^^M p4Mzlu dz 

£M^ M ^, (z)udz 

- 1 p=0 

= 1 - \A\ +[ L x , u (d, (3' (z),f x (z)s) dz 
J A 



(29) 



where \ A\ is the measure of A and L x ^ u (d, (3, s) is the Laplace transform of f x ^ u (d, (3, z). By using 
the properties of the Laplace transform and by taking its inverse, we finally get 

f x , x (d,f3,z) = (l-\A\)5(z) + [ -J—f X)U (d,-^- -^) dz. (30) 

We can rewrite the second term of (1301 by defining the cumulative density function 

G x (y) = jL\ {zeA \f x (z)<y}\ 

for y > 0. By using the corresponding probability density function, g x {y), and Lebesgue integration, 
we can rewrite (1301 as in (fTTT) . 

Appendix C 
Corollary 15. II 

From the result in (ITTb and from the assumption / x (z) > Vz £ % (i.e., |.A| = 1), we have 

h, x (d, p,z)= / /a,« d, -, - dy 

Jo 2/ \ y yj 

Then, from the definition of f x >(z) given in (fT2l it follows that G x '{y) = G x {cy) and, by consequence, 
9x'(y) = cg x {cy). Therefore, from (fT3l we have: 

fx,Ad,P,z) = (l-c)d(z) + c r^Mf^fd,^-,-) dy 



o y \ y y 

h W ^ 2 f°° 9x(cy) f f,PA A 

Jo y \ y yj 

(l-c)5(z) + c 2 f x>x (d,cf3,cz) 
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Appendix D 
Corollary [53] 

From the expression of the moments given in Theorem l5.1l and the results in (9), it is easy to show 
that for uniformly distributed phases, we have: 



lim M p d j3 x = 1 



for any p > 0. It immediately follows that 



lim f\, u (d,/3,z) = 5(z - 1) 

/3->0 



where S(z) is the Dirac's delta function. By applying this result to (fTTb . we get 

^ hJdAA = (l-\A\)5(z) + \A\ r?Ms(^-l) dy 

/3^o v v vJ Jo y \y ) 



(i - |^|) 5(z) + \a\ / ^ z,w) 5{w - 1) dw 



o w 

(l-\A\)5(z) + \A\g x (z). (31) 
Appendix E 



Corollary 15 .21 

By using the definition of the 77-transform and the result in (fTTT) . we obtain: 

i] x {d,M) = ^[(tA + i)" 1 

c 00 1 

— -f\,x{d,P,z)dz 

7^ + 1 



7^ + 1 7 y Jo 72 + 1 V y yy 

1-1-41 + 1-41 9x{y) I" l -—fx,u(dA,^) ^dy 

Jo Jo 7F + 1 V y ) 

\-\A\ + \A\j g x {y)vu(dA,iy 



dy. (32) 

Then, by considering that MSEqo = rj x (d, P,j/j3), the expression of the asymptotic MSE immediately 
follows. 

Appendix F 
Corollary [53] 

In Appendix 15.31 we have shown that lim^o f\, u (d, (3, z) = S(z — 1). From the definition of the 
77-transform, it follows that lim^o r] u (d, (3, 7) = and lim^o Vu(d, f3,j/(3) = 0. Thus, from ([141 ) 
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we have: 



lim rj x (d,P, 7) = 1 - \A\ + \A\ [ g x (y) ( lim i] u ( d, -,^y ) ) dy 
0-+° Jo W° V V J J 



00 



l-\A\ + \A\ I -^Idy 



7?/ + 1 

= 1-1^1 + 1^1^(7) (33) 
where we defined ^(7) = / °° *^pf dy. As a consequence, 

\hn Vx (d,P, 7 /P) = 1-\A\ + \A\ T g x (y)(hm Vu (d^,^f)) dy 

Jo \p^o v y p J J 

= l-\A\. (34) 
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